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We present a comprehensive analysis of high resolution neutron scattering data involving Neutron 
Spin Echo spectroscopy and Spherical Polarimetry which confirm the first order nature of the helical 
transition and reveal the existence of a new spin liquid skyrmion phase. Similar to the blue phases of 
liquid crystals this phase appears in a very narrow temperature range between the low temperature 
helical and the high temperature paramagnetic phases. 



I. INTRODUCTION 

MnSi is one of the most investigated chiral magnets 
and, at least in theory, a very simple realization of chiral 
magnetism. The Ginzburg-Landau Hamiltonian contains 
three hierarchically ordered interaction terms with well 
separated energy scaled. It is therefore possible to distin- 
guish between the different contributions to the ground 
state. The strongest ferromagnetic exchange interaction 
term fixes the spins at the longest range. The weaker 
Dzyaloshinskii-Moriya (DM) term arises from the lack of 
inversion symmetry of the B20 lattice structure and ro- 
tates the spins at the intermediate scales. And the weak- 
est anisotropy term fixes the directions of the spins on the 
crystallographic lattice. The magnetically ordered state 
below Tc 29 K is a single domain left handed helix 
with a period of '^175 A . All magnetic moments are per- 
pendicular to the propagation vector r that points along 
the <111> crystallographic directions'^. In the helical 
phase small angle neutron scattering shows well defined 
Bragg peaks of the same intensity at all equivalent mag- 
netic reflections |tiii|, with r= |riii| =0.036 A~-^. 

MnSi is a weak itinerant magnet with an ordered mag- 
netic moment of only 0.4 /i^, a fraction of the effective 
magnetic moment of 1.4 jub determined in the paramag- 
netic phase^. The strong magnetic fluctuations, which 
are due to the weak itinerant magnetism and the vicinity 
of a magnetic instability, exist not only above Tc^ but 
also persist in the low temperature ordered phase^. The 
magnetic fluctuations lead to an enhanced effective elec- 
tron mass, a broad specific heat maximum reminiscent 
of the specific heat of spin liquids, frustrated magnets or 
spin glasses^ and broad features on thermal expansion 
and ultrasound measurements, which almost completely 
mask the helical transitionF. 

One intriguing feature of MnSi is that magnetic corre- 
lations above Tc appear not only around the positions 
in reciprocal space of the helical order but spread homo- 
geneously over the whole surface of a sphere with radius 
r emerging as a powder diffraction like ring on the two- 
dimensional small angle neutron scattering spectra^. If 



the neutron beam is polarized the rings reduce to half- 
moons due to the interaction between polarized neutrons 
with the helical correlations as explained below. 

This unconventional feature occurs in a limited tem- 
perature range above Tc and is reminiscent of scattering 
patterns from cholesteric liquid crystals suggesting the 
existence of similar textures also in magnets. However, 
the chiral molecules of liquid crystals are rods whereas 
magnetic moments are vectors. This additional topolog- 
ical constraint leads to the formation of domain walls 
rendering tt disclinations, i.e. the low order line defects 
commonly found in liquid crystals, energetically unfa- 
vorable in magnetJ^. On the other hand, domain wall 
formation is minimized by higher order skyrmion-like 27r 
disclinations. 

Skyrmions were introduced in the early 60s by 
Skyrme^^ to bridge the gap between waves and parti- 
cles in the particle-wave duality description. The exis- 
tence of these soliton-like quasiparticles was alleged in 
semiconductors under high magnetic fields^' and their 
topology corresponds to that of the blue phases of liq- 
uid crystals^^. Periodic skyrmion lattices have recently 
been identified under a magnetic field in the A-phase of 
MnSP^ or directly seen by Lorentz transmission electron 
microscopy in Feo.sCoo.sSi^^ or FeGe^^. 

It has been suggested that the unconventional diffuse 
scattering above Tc reflects the existence of skyrmion- 
like textures, which would form spontaneously without 
any structural defects or an external magnetic field^ sta- 
bilized by spin stiffness or higher order termJ^^^. This 
novel ground state of condensed matter would result from 
the isotropically twisted spin liquid phase and as in liquid 
crystals it would occur at a very restricted temperature 
range just above the long range ordered helical phase. In 
the following we present high resolution polarized neu- 
tron scattering experimental results, which strongly sup- 
port this hypothesis. 
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FIG. 1: Principle of the Polarimetric Neutron Spin Eclio teclinique. Tlie neutron wavelengtli information is encoded in the 
first precession arm between two 7r/2 flippers. The beam is re-polarised before entering the zero-field region of Cryopad. The 
incident and scattered polarisation vectors are set by Cryopad. After scattering, the second precession arm, between two 
additional 7r/2 flippers, encodes again the wavelength and the echo is measured after the analyser. 
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FIG. 2: Temperature dependence of the logarithm of the in- 
tensity at the position of the helical peak (riii= 0.036 A~^). 
A intensity jump of almost one order of magnitude defines 
the helical transition Te- 



ll. EXPERIMENTAL METHOD 

MnSi crystallizes in the cubic P2i3 (T4) structure with 
the Mn atoms occupying the 4a site for x = 0.138. The 
lack of inversion symmetry in the crystallographic struc- 
ture is at the origin of the DM interaction, the chiral 
magnetism and the helical magnetic structur e. T he ex- 
periments were done on a well characterized^^ single 
crystalline sample with a thickness of 2 mm and a diam- 
eter of 20 mm, cut from a large single crystal grown at 
Ames Laboratory. The lattice constant was 4.558 A. The 
sample was oriented with a <110> direction vertical, so 
that four 111 reflections and two 110 reflections were ac- 
cessible in the horizontal scattering plane. 

The polarized neutron scattering experiments were 
carried out with the Neutron Spin Echo spectrometers 
INll and IN15 at the Institut Laue Langevin, Greno- 
ble, France. On IN15 the wavelength was 8 A and on 



INll 6.5 A. At both instruments the incoming neutron 
beam had a monochromatization of 15% FWHM. The 
XY position sensitive detectors of both spectrometers 
covered an angular range of ~3x3 deg and the resolu- 
tion in momentum transfer Q was ~ 0.005 A-i FWHM. 
On INll the Q dependence of the relaxation and the 
intensity was analyzed by taking slices of the detec- 
tor at constant distance from the 000 point and also 
from the surface of the sphere with radius r leading to 
-0.005 A"^ < q = \Q\-r < 0.018 A"\ The measure- 
ments were carried out at one of the four equivalent mag- 
netic reflections OOO+rm. On INll it was also possible 
to measure at a OOO+rno point, where the correlations 
do not develop to a Bragg peak. 

Fig. |2] shows on a log-lin scale the temperature depen- 
dence of the neutron intensity at rm. The helimagnetic 
phase transition is marked by an intensity jump of more 
than one order of magnitude within less than 0.1 deg. 
The data above and below Tc follow power laws similar 
to those found previousl}El and which were at the origin 
of a long standing controversy on the order of the tran- 
sition. These power laws, however, do not account for 
the intensity jump, which confirms the first order phase 
transition seen by specific heat(^. 

The intensity jump defines very accurately the transi- 
tion temperature is Tc = 29.05 ± 0.05 K. These data 
were collected on INI 5. A different cryostat and ther- 
mometer were used on INll and the jump was observed 
at Tc = 28.6 ±0.05 K. In order to compare the data sets 
the results will be plotted against T-Tc or the reduced 
temperature e =(T-Tc)/Tc. 

The helical Bragg peaks are the fingerprint of the he- 
lical phase. They are elastic (energy transfer ficj = 0, in 
the following we will often use oo to also designate the en- 
ergy transfer without explicitly using h) and has a Gaus- 
sian lineshape. In contrast, the fluctuating paramagnetic 
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phase above Tc has finite correlations and the scattering 
function S{Q,uj) is a superposition of Lorentzians^ ^. In 
the quasi-elastic limit, where the energy transfer is much 
smaller than the energies of the incoming beam and of 
the sample, the neutron scattering cross section becomes: 



dQ dE 



(7 C 

oc S{Q,uj) oc ^ 



.2 r2. 



(1) 



with C the Curie constant, = 1/^, where ^ is the char- 
acteristic correlation length, and T the HWHM energy 
linewidth. S{Q^uo) is therefore the product of the static 
structure factor, which has the Ornstein-Zernike (OZ) 
form : 



S{Q) = [ S{Q,uj) duj oc - 
and of the dynamic structure factor: 
syQ^uj) — (X 



C 



S{Q) 



Y2 



(2) 



(3) 



The dynamic measurements were performed by Neu- 
tron Spin Echo spectroscopy, which uses the Larmor pre- 
cession of the neutron spin in a magnetic field as a clock 
to measure with very high accuracy the difference in neu- 
tron velocities before and after the scattering process at 
the sample. The changes in the neutron energy due to 
inelastic scattering affect the neutron velocity and conse- 
quently the amplitude of the Larmor precessions. In this 
way, the energy transfer is measured directly by circum- 
venting the intensity resolution limitations of the Liou- 
ville theorem. For this reason NSE reaches very high res- 
olutions while maintaining the high intensity advantage 
of a beam that is only 10-20% monochromatic. The high- 
est energy resolution in neutron scattering is presently 
reached by IN15, which accesses energies as low as some 
neV corresponding to motions with characteristic times 
reaching the /is. 

At the quasi-elastic limit, which is valid in most NSE 
experiments, the amplitude of the NSE Larmor preces- 
sions is directly proportional to the intermediate scat- 
tering function /(Q,^), the Fourier transformation of 
s{Q,uj^P^. The comparison between NSE and inelastic 
neutron scattering is therefore straightforward through a 
Fourier transformation. 



IiQ,t) = 



f S{Q,uj) cos{ujt) dio _ ?ft{S{Q,t)) 



J S{Q,oj) du) 



S{Q) 



(4) 



Consequently, the Lorentz function of eq. [3] Fourier 
transforms to the exponential NSE decay exp(— t/to) 
with to = 1/r leading to t^lns] = 0.658/r[/ieV]. With 
this relation it will be possible to compare our NSE re- 
sults with the triple axis neutron spectroscopy (TAS) 
measurements from the literature^ 

This so-called Larmor labeling requires polarized neu- 
trons and some polarization analysis features are integral 




FIG. 3: Temperature dependence of the chiral fraction deter- 
mined by neutron polarimetry as explained in the text. 



part of NSE^^. The pulse sequence, 7r/2 fiip-precession- 
TT fiip (at the sample)-precession-7r/2 fiip, is that of the 
classical Hahn sequence of NMR SE. "Normal" param- 
agnetic scattering acts as a tt flipper and gives an echo 
without the otherwise obligatory tt flipper, which leads 
to a straightforward and unambiguous separation of the 
magnetic and nuclear signals. On the other hand, due 
to Larmor precessions, the neutron beam is depolarized 
at the sample with all neutron spins evenly distributed 
in the precession plane. For this reason, in the pres- 
ence of chiral and/or nuclear-magnetic interference terms 
the analysis of the experimental results is complex. The 
way out is the Polarimetric Neutron Spin Echo set-up, a 
variant of Intensity Modulated NSE^^, which combines 
the precession fleld areas required for Neutron Spin Echo 
spectroscopy with CryopacP^. As shown schematically 
by flg. [l] the precessions are stopped before the sample 
by an additional 7r/2 flipper. The neutron beam is then 
re-polarized with a compact polarizer to allow full analy- 
sis of the scattered polarization vector with the zero fleld 
chamber of Cryopad. The precessions resume at the sec- 
ond branch of the spectrometer after a second additional 
7r/2 flipper and the echo signal is recovered at the neu- 
tron detector after the analyzer. 

This set-up is now implemented at IN 15 and can be 
used to observe chiral fluctuations with unprecedented 
accuracy and resolution both in energy (time) and mo- 
mentum transfer (space). The Polarimetric NSE and 
spherical polarimetry measurements were performed with 
a third generation Cryopad, a sensitive zero fleld po- 
larimeter which controls the polarization of the incoming 
and the scattered beams^^^^. A combination of /i- metal 
and Meissner shields reduces the residual magnetic fleld 
at the sample position down to ~0.1 /iT and controls the 
direction of the polarization vectors with an accuracy of 
better than 1 deg. The beam polarization was 96% cor- 
responding to a flipping ratio of 45. 



4 



III. POLARIZED NEUTRON FORMALISM - 
SPHERICAL POLARIMETRY 

Before presenting in detail the experimental results we 
will introduce the interaction between a polarized neu- 
tron beam and the magnetic helix of MnSi following the 
formalism developed almost simultaneously by Blume^^ 
and Maleye\^^ in 1963. The helix is described by two 
orthogonal vectors si and 52, which have the same am- 
plitude and are perpendicular to the helix propagation 
vector r. The magnetic structure factor for a reflection 
K is given by: 



M{K) = ^^rhj f{K)exp{2mK ■ 



(5) 



with rhj 



■ jii{si — is2) exp{2m r ■ fj) exp{i(j)j) 



where ja is the amplitude of the Mn moments, f{K) the 
magnetic form factor of the Mn atoms, fj an atomic po- 
sition, and an arbitrary phase angle. As the mag- 
netic interaction vector M±{K) = K x M{K) x K is the 
projection of the magnetic structure factor onto a plane 
perpendicular to the scattering vector Q = i?, it is con- 
venient to choose the set of orthogonal polarization axes 
with: 

HQ, 

• z perpendicular to the scattering plane and 

• y completing the right-handed cartesian set. 

In the geometry of the present experiment y=n with h 
the propagation vector of the incoming neutron beam. 
For the polarimetric measurements the MnSi sample was 
oriented so that Q = fm and z\\ (1, 1, 0). For this reflec- 
tion the intensity of the diffracted beam was proportional 
to: 



a = M^Ml + P' ■ ^(Mx X Ml 



(6) 



where is the polarization of the incoming beam. The 
second term is the chiral part characteristic of the spiral 
structure and it adds to the conventional first term for 
P' = Q i.e. when the incident polarization is parallel 
to Q. On the other hand, the reflection is completely 
extinct for = —Q. Following Blume's equations^^, the 
scattered polarization for MnSi is given by: 

Per = -P'{M^Ml) + 23?((P • Ml)M^) - ^(M^ x Ml) 

(J) 

The first two terms form the 'trivial' magnetic part lead- 
ing to P = —Q {Q ■ P') for isotropically distributed elec- 
tronic magnetic moment^^. The chiral third term cre- 
ates a polarization antiparallel to Q independently from 
the polarization state of the incoming beam. Both equa- 
tions [6] and [7| lead to the rigorous and accurate determi- 
nation of the chiral component ^(M^ x M^) by polarized 
neutrons. The more general form of eq. [7|is: 



with P the polarization transfer tensor and Pt = 
— ^{M±_ X Ml) the polarization created by the chiral 
sample. It is then useful to have the incident polariza- 
tion along —X, X, y and z and to measure the outgoing 
polarization components along ±x, ±^ and =bz for each 
configuration of the incident polarization respect ivel}i^. 
This measurement procedure is characteristic of spherical 
neutron polarimetry (SNP) and determines the polariza- 
tion matrix : 



P' . 



p. . p' 



PI 



with {ij) e {x,y,z} (9) 



where the denominator |P| corrects for the finite effi- 
ciency of the neutron polarizer-analyzer system and the 
intrinsic imperfections of the experimental setup. For nu- 
clear (coherent) scattering ¥ij = Sij, i.e. 1 for z = j and 
otherwise: 



- nuclear 



1 
1 
1 



For an ideal paramagnet all matrix elements are zero ex- 
cept = — 1: 



L para 



-10 





In MnSi according to eq.[7|the scattered beam will always 
have a polarization antiparallel to the helix propagation 
vector, i.e. antiparallel to x and Q. For this reason the 
first row of the matrix will be non-zero, ideally ¥i^x = — 1, 
'^i,yorz = with (i) G {x, z}. In its most general form 
the chiral matrix can be written as: 



^ch 



iral 



-1 r]C r]C 





where ( determines the chirality of the helix: = +1 for 
right and = — 1 for left handed chirality. rj measures the 
fraction of the dominant chiral domain: 77 = 1 for a single 
domain and 7^ = for equally populated chiral domains 
or for the disordered paramagnetic state, in which case 

^chiral reduces to Ppara- 

The polarization matrix of MnSi below Tc is that of 
an ideal chiral left handed single domain structure. At 
Tc-4 K we found: 



-4K = 

-1.000 ±0.001 
-0.007 ±0.001 
0.054 ± 0.003 



-0.995 ±0.001 
0.016 ±0.002 
0.055 ±0.002 



-1.002 ±0.001 
-0.007 ±0.002 
0.062 ± 0.002 



P = PP' ^ Pt 



(8) 



Just above Tc, the intensity at rm drops dramatically, 
which leads to higher counting times and error bars. Nev- 
ertheless, the matrix remains unaffected and is that of a 
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perfect left handed single chiral domain up to ~Tc+l, 
e.g. at Tc+0.4 K: 



+0.4K 



-1.03 ±0.03 
-0.01 ±0.05 
0.07 ±0.05 



-0.99 ±0.03 
-0.00 ±0.05 
0.04 ±0.01 



-0.98 ±0.03 
-0.03 ±0.05 
0.07 ±0.05 



Above ~Tc±l the non-diagonal elements decrease slowly 
and chirality remains finite even well above Tc- E.g. at 
Tc±4 K we found: 



-1.02 ±0.02 
-0.01 ±0.1 
0.02 ±0.1 



-0.4 ±0.1 
-0.08 ±0.1 
0.07 ±0.1 



-0.3 ±0.1 
-0.05 ±0.1 
-0.08 ±0.1 



All data have been corrected for the background deter- 
mined from polarimetric measurements at 50 K. Fig. |3] 
shows that F^y and F^z are within the error bars unaf- 
fected by Tc and start to decrease only above ~Tc±l K. 
This result is confirmed by the intensity of the reflection 
when P' \\ ± X, in which case eq. [g] becomes: 



-^para ± -^chiral 



with 



(10) 



para 



M^Ml and N, 



chiral 



where the sign in front of the chiral contribution, NchiraU 
depends on the direction of P' (parallel or antiparallel) 
with respect to Q. The intensities in the spin flip channel 



P' II Q 



^^bck 



x,x — ■'-^ pn\_ 



SF 



bck 



' -^chiral 
■^chiral 



(11) 



The chiral fraction N chiral / {N chiral ^ N para) can then be 
deduced from the ratio {Nx-x—N-x^x)/{Nx-x^N-x,x — 
2*A/'^^), with the background determined at 50 K. 
This method minimizes the error bars and systematic er- 
rors because both N^-x and N-x,x are measured under 
identical conditions. The deduced chiral fraction is also 
shown in flg. [3] and coincides with the non-diagonal ma- 
trix elements. 

Thus MnSi remains completely single domain chiral 
up to TJ^ ^Tc±l K, i.e. in the complete temperature 
range, where the isotropic scattering at the surface of the 
sphere with radius r is observed. We identify TJ^ as the 
temperature, where the theoretically predicted skyrmion 
spin liquid phase sets in, and the analysis of the dynamic 
and static correlations in the following sections will con- 
flrm this hypothesis. 



IV. NEUTRON SPIN ECHO 

The dynamics of MnSi was thoroughly investigated 
by triple axis spectroscopy (TAS) in the mid-80&^. It 
was shown that strong magnetic correlations exist up 
to room temperature and that the magnetic excitation 
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FIG. 4: Dynamic correlations at the position of one of the 
helical Bragg peaks rm measured by polarimetric NSE for 
the Fz,x term of the polarization matrix (eq. [oj. 
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FIG. 5: Dynamic correlations at rm measured by classical 
NSE. 



spectrum is well described by the Moriya-Kawabata the- 
ory for weak itinerant ferromagnets. The TAS measure- 
ments, however, did not have the resolution in energy 
transfer required to analyse the magnetic fluctuations 
close to Tc. With Neutron Spin Echo spectroscopy we 
can reach the required resolution, follow very accurately 
the slowing down of the magnetic fluctuations and com- 
plement the TAS experiments. The fact that MnSi is 
completely chiral above Tq was an unexpected outcome 
of the polarimetric measurements and also a posteriori 
the justiflcation for combining NSE and Cryopad to po- 
larimetric NSE. 

With the Cryopad on IN 15 we were able to measure 
separately the relaxation of the diagonal ¥xx and of the 
crossed chiral terms ¥yx and F^x- The corresponding in- 
termediate scattering functions I(q,t) were obtained by 



6 




Tc+1.4K 



Tc+0.6 K 



0.001 0.01 



Fourier time [ns] 

FIG. 6: Dynamic correlations at rno, where the spectra are 
not contaminated by the hehcal Bragg peaks. 



normalizing the spectra against the resolution measured 
below Tc, typically at 25 K. Fig. [4] shows polarimet- 
ric NSE spectra for the crossed term ¥yx at rm. This 
was the first time that magnetic chiral fluctuations were 
measured directly and independently from the trivial di- 
agonal if'xx) part. 

The polarimetric NSE spectra display a purely expo- 
nential relaxation superimposed on an elastic term: 



= exp(— t/to) + (1 — a) 



(12) 



As seen in fig. |4]the elastic fraction (1 — a) evolves from 
^15% to 100% within 0.15 K following the fast increase of 
the intensity displayed in fig. |2] and masking the dynam- 
ics at Tc. This behaviour is confirmed by standard NSE 
measurements (fig.[5|, which however do not extrapolate 
to 1 at t^O. The comparison with polarimetric NSE 
identifies this as an artefact due to the chiral magnetic 
scattering and not to additional dynamic components. 

At Tiio the quasi-elastic scattering does not develop 
to a Bragg peak and the NSE spectra were not contami- 
nated by an elastic contribution. Consequently the fluc- 
tuations could be followed even below Tc (fig. [6|. The 
decay is exponential with the same to and deduced F as 
at Till. 

All relaxation times and linewidths at = are plot- 
ted against temperature in fig. [7| The diagonal {xx) and 
crossed terms {yx and zx) measured in the polarimet- 
ric NSE mode are also included. All data, even those in 
the closest vicinity of Tc, fall on the same curve. Con- 
sequently, the relevant parameter for the fluctuations is 
the distance from the sphere with radius r, not the Bragg 
peaks Till of the helical phase. 

We will now discuss the g-dependence of the relaxation. 
Close to Tc the linewidths take their lowest values at 
the surface of the sphere with radius r. We found that 
F increases both ioi \Q\ > r and |Q| < r and the low-q 
points in fig. [s] are the average for q=±0.03 A How- 
ever, above ^Tc+1.4 K the q-dependence of F flattens 
and it is no longer possible to identify the position of the 
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FIG. 7: Temperature dependence of the characteristic times 
to (a) and the deduced Lorentz linewidths (b) of the fluctu- 
ations at Till and tho measured by classical (paramagnetic) 
and polarimetric NSE. The lines are guides to the eye. 
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FIG. 8: q dependence of the dynamic linewidth F. The open 
and closed symbols represent the values at rm and rno re- 
spectively. At the highest temp eratu re {r^ 3 1 K) there is ex- 
cellent agreement with literatur^^^ and eq. 17b (red dotted 
line) . 
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FIG. 9: Magnetic signal, proportional to S{q), around 
the Till. The continuous lines are the best fits: simple 
Gauss function at Tc, superposition of a fluctuating Lorentz 
(Orstein-Zernike) and an elastic Gaussian part at Tc < 
T < Tc + 0.2 K and a fluctuating Lorentz function above 
Tc +0.2 K. 



minimum. 

Close to Tc T{q = 0) levels off at - 0.64 /ieV, which 
implies that the associated correlation length does not 
diverge at Tc- The results will be discussed in the frame 
of dynamic scaling after the determination of the corre- 
lation length in the following section. 



V. CORRELATION LENGTH 

The Q dependence of the magnetic static structure 
factor S{Q) was analysed in the configuration N^^-x of 
eq. [TT] where the magnetic intensity is maximum and the 
background correction negligible for T<30 K. 

Fig. [9] shows the scattered neutron intensity, which 
is proportional to S{Q)^ around rm in a log-lin scale. 
At Tc the points are best fitted with by a Gauss 
function, with = 6.1110-^ ±310-^A-^ leading to 
1.4410-2 ± 10-"^ FWHM, i.e. significantly broader 
than the Q-resolution of the instrument. Well above Tc, 
the data are well described by the OZ function of eq. |2] 
convoluted with the Gauss resolution function. 

Just above Tc, however, neither the Gauss nor the 
OZ functions fit satisfactorily the experimental data. In- 
stead, the best fit is obtained by a superposition of a 
fluctuating OZ and an elastic Gauss with relative weights 
fixed from the NSE spectra, as shown by fig. [lO] We note 
that all three fits in fig. 10 involve the same number (two) 
of independent parameters: the total intensity and the a 
for the Gauss or k for the Lorentz and Gauss+Lorentz 
fits respectively (in this case the width of the Gauss was 
fixed to that found below Tc). 



FIG. 10: Comparison of the quality of fits for a simple a simple 
Lorentz (black dotted line) , Gauss (blue dashed line) , and the 
superposition of a Lorentz and a Gauss (red continuous line) 
at Tc+0.05 K and rm 
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FIG. 11: Plot of the correlation lengths (^m and ^110) mea- 
sured at Till and rno as a function of temperature. 



Similarly to the NSE spectra, the data analysis is eas- 
ier at Tiio, where there is no contamination from Bragg 
peaks. At this position of the reciprocal space a simple 
Lorentzian describes the scattered neutron intensity at 
all temperatures, even below Tc- All correlation lengths 



are plotted varsus T-Tc in fig. 11 and the deduced values 
of Hi are in excellent agreement with previous published 
data^ (fig.[l2|. 

At Tiio the correlation length ^no levels off at about 
100 A , i.e. about half the pitch of the helix. At rm the 
correlation length increases considerably close to Tc 
following a power law of the reduced temperature e sim- 
ilar to those found at second order phase transitions and 
illustrated by fig. 13 The continuous line in the figure 
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T^^Q after Grigoriev et al. 
T-,-,-, after Grigoriev et al. 



T-Tc [K] 



2.5 



FIG. 12: Plot of n measured at rm and rno- The closed 
symbols are the present work. The open symbols are from 
Grigoriev et al.^. 
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FIG. 13: Plot of the correlation legth determined at rm 
versus the reduced temperature e on a log-log scale. The 
continuous line is the best power law fit ^ = 12 e-°-^ A. The 
black dotted line shows the extrapolated values from the high 
temperature TAS work of Ishikawa et al.^. 



corresponds to: 

^111 = a e"^ with a = 12 ±2 A and i/ = 0.5±0.2 (13) 

The black dotted line represents ^ = 5.6 e~^"^ A, the ex- 
trapolated curve from the high temperature TAS data of 
Ishikawa et al.^. It is remarkable that this extrapolation 
from very high temperatures (100 K - 300 K ) is only a 
factor 2 off our results. 

VI. DYNAMIC SCALING 

Dynamic scaling relates T to z^: (or inversely to to ^) 
through a homogenous function: 



(14) 



with the dimensionless ratio n/q defining the critical 
(/^/(7 <C 1) and hydrodynamic {n/q ^ 1) regimes 
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FIG. 14: Temperature dependence of the ratio r(q=0)//€iio 
revealing the existence a broad flat minimum, where r(q=0) 
oc A€iio between Tc and T'q. 
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Tr+Z.3 K 




FIG. 15: Scaling plot of the reduced dimensionless quanti- 
ties r/r(Tc) versus Kno/q. The black dotted line shows the 
extrapolat ed va lues from the previous high temperature TAS 
data (eq. 17b). The red dashed line is the Resibois-Piette 



function for antiferromagnets and the green continuous line is 
a guide to the eyes through the experimental data for Tc <T 



respectivel}^^. At the critical limit, which is always 
reached at k,=0, the linewidth and relaxation times re- 
flect the volume probed by the measurement : 



to(<7,e 
r(g,e 



0) (X q- 
0) (X q^ 



and 



(15a) 
(15b) 



with z the dynamic exponent. For (7 = Q as in ferromag- 
nets z = 5/2, whereas for antiferromagnets and MnSi, 
where the relevant parameter is q = \Q — z = 3/2. 

At the hydrodynamic limit, q = and > 0, the dy- 
namics is governed by the correlated volumes through the 
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power laws: 



to{q = 0) cx and 
T{q = 0) cx K.^ 



(16a) 
(16b) 



The previous extensive study of spin fluctuations in 
MnSi with TAS^ showed that the linewidth is equally 
well described by two universal scaling functions over an 
extremely large temperature scale, from ~30 K up to the 
room temperature: 



(17a) 
(17b) 



with = 0.0325 e 1"^ 

Ao=19.6 mevA^/^. Equation 
the Moriya-Kawabata theor;^ 



17a 



=50 meVA^ and 
is derived from 
for weak itinerant 
ferromagnet^^ whereas eq. |l7b| has the form expected 
for dynamic scaling. The red dotted line going through 
the data at Tc+2.3 K in fig [S] corresponds to the calcu- 
lated values from eq. |17b| and is in excellent agreement 
with our experimental results. 

At high temperatures, the momentum transfer Q 
is the relevant parameter and z = 5/2. Close to Tc, 
however, magnetic correlations build up at r and the 
relevant parameter crosses over to q = |Q — r|, which 
is used in this work. Consequently also the dynamic 
critical exponent should cross over from z = 5/2 to 
z = 3/2. 

If the transition were of second order, an analysis 
along the lines of dynamic scaling would imply critical 
slowing down. In this case Tiii{q = 0) should decrease 
continuously to zero at Tq following the divergence 
of ^111. In contrast, Tiii{q — 0) remains finite and 
strictly the same as Vw^ici = 0). Moreover, fig. [t] and 
fig. [12] show that close to Tc both Fin and Fno are 
roughly proportional to a^hq, not to the "critical" hzm. 
Consequently, we do not observe the critical slowing 
down and rapidly changing dynamic behaviour expected 
for second order phase transitions. This is underlined by 



fig. [14] where the ratio between all F and the calculated 
Hiiio levels off to a broad plateau between Tc and TJ^, 
pointing towards z ^ 1. 



Fig. [15] shows the interdependence of the reduced 
dimensionless quantities F/F(Tc) and /^no/q- Following 
the scaling assumption of eq. [T4| all experimental data 
should collapse on a universal curve. This is obviously 
not the case. Close to Tc both F and tzno do not 
vary with temperature and for this reason the data 
accumulate on the green continuous line of fig. [Tsj 
which is significantly different from the Resibois-Piette 
function for antiferromagnetP^. Above ~TJ^ there is a 



n J? 



Paramagnetic Phase 
partially chiral 




FIG. 16: Temperature dependence of (a) the intensity at the 
helical Bragg peak rm (b) the chiral fraction, (c) the ratio 
T/kiio and (d) the specific heat from Stishov et al.^. These 
experimental results confirm the existence of a new completely 
chiral but fluctuating skyrmion liquid phase between Tc and 
T'c- 



fast change of behaviour crossing over to eq 17 and the 



black dotted line in fig. 15 
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VII. DISCUSSION 

These high resolution neutron scattering data con- 
firm the first order nature of the helical transition in 
MnSi seen b y spe cific heat and ultrasonic attenuation 
measurements^!^. In favour of the first order phase tran- 
sition are: 

a. the sharpness of the transition (fig.[2|, even though 
the helical Bragg peaks are not resolution limited. 

b. the coexistence of high and low temperature phases 
between Tc and Tc+0.2 K. At rm the NSE spectra are 
the superposition of a fluctuating and elastic part (fig. |4] 
and [5| and the static structure factor is best described 
by the superposition of a Lorentzian (fluctuating) and 
Gaussian (elastic) contributions (fig. [lO]). 

c. the absence of critical slowing down at T^. The 
slow change of the dynamics scales with the non diverg- 
ing correlation length ^hq. 

Fluctuations are present above Tc and can indeed 
be approximated by a power law (eq. 13 and fig. 13). 



For this reason, a selection of experimental data can al- 
ways be bent to fit into the second order phase transition 
scheme by introducing crossovers, as it was recently done 
by Grigoriev at al.^^. However, second order phase tran- 
sitions are global transformations involving all parame- 
ters of the system, and this is not the case in MnSi. 

We will now proceed to the discussion of the fluctu- 
ating phase between Tc and TJ^. Fig. 16 recapitulates 



the most important findings of this work including the 
specific heat measured by Stishov et al.^^. The first or- 
der helical transition Tc is marked by the sharp peak in 
the specific heat (fig.[l6]i) and the Bragg peak (fig. 16 1). 
However, the helical transition has no effect on 77, the de- 
gree of single domain chirality. Comparison with theory 
shows that this is not a trivial effect, r] can indeed be 
calculated in the frame of a simple mean field model^ : 



V = 2QT/{Q^ + T^ + l/e) 



(18) 



The measured lead to the dotted green line of 

fig. [16]3, which is significantly different from the experi- 
mental data. Consequently, the simple assumption of an 
unpinned and fluctuating helical phase is not confirmed 
by the experiment. In contrast, 77= 1 implies a high de- 
gree of chiral short range order for the fluctuating phase 
between Tc and Tj^. 

TJ^ is exactly the theoretically predicted temperature 
for the formation of a skyrmionic ground state in MnSi 



Fig. 16 shows that this is not a singular event but part of 
a global transformation involving the dynamics as well as 



a broad maximum in the specific heatP^I. Besides specific 
heat, thermal expansion, sound velocity, sound absorp- 
tion as well as resistivity measurements point towards 
a non trivial spin liquid phase^. Skyrmions account for 
the inherent topology of this completely chiral spin liq- 
uid phase, which gives a scattering with rotational sym- 
metry. Skyrmionic textures are energetically favourable 
at short distances, which explains the slow decrease of 77 
and its finite value well above Tc- The correlation length 
^^110 gives a natural measurement of the skyrmion radius, 
which peaks are about half of the helix pitch close to Tc, 
in agreement with the skyrmion hypothesis. 

The Bragg peak of fig. [16^ is a direct measurement 
of the helical order parameter. We speculate that r] in 
fig. ^Ejp reflects the order parameter of the skyrmion spin 
liquid phase. The topology of skyrmions would explain 
the particular dynamics and the low dynamic exponent 
of z ~ 1. This skyrmion liquid phase would then be 
the magnetic equivalent of the blue phases in liquid crys- 
tals, which occur at a very limited temperature range 
above the helical phase. In this way the phase diagram 
of MnSi becomes very similar to that of cholesteric liquid 
crystals, emphasizing the parallel between magnetic and 
structural ground states. Similarly to liquid crystals, the 
transition to the helical state is of first order due to the 
different topologies of the helical and skyrmionic phases. 

In summary, we combined high resolution neutron spin 
echo spectroscopy and spherical polarimetry to obtain a 
consistent picture of the spin liquid phase of MnSi above 
Tc- The results evidence a first order transition between 
the helical phase and a skyrmion liquid phase. This com- 
pletely (single domain) chiral and strongly fluctuating 
new state of matter occurs in a narrow temperature range 
above the helical phase of MnSi. 
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